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ABSTRACT 

The ultraviolet (UV) spectra of 16 blazars (< z >~ 1) from the archives of the Hubble 
Space Telescope Faint Object Spectrograph have been analyzed in order to study in a sys- 
tematic way the properties of their broad UV emission lines. We find that the luminosities of 
the most prominent and intense lines, Lya and C IV A 1549, are similar to those of normal 
radio-loud quasars at comparable redshifts. However, the equivalent widths of blazar lines are 
significantly smaller than those of radio-loud quasars. Therefore, while the intrinsic broad line 
region luminosity of blazars appears to be indistinguishable from that of radio-loud quasars, 
their continuum must be comparatively higher, most probably due to relativistic beaming. We 
have combined the UV luminosities of the de-beamed continuum with the emitting gas veloc- 
ity to derive estimates of the masses of the central supermassive black holes. The size of the 
broad line region was computed in two ways: 1) via an empirical relationship between UV 
continuum luminosity and broad line region size, and 2) through the external photon density 
required by blazar models to reproduce the inverse Compton components observed at gamma- 
rays. The second method yields significantly different results from the first method, suggesting 
that it provides only a very rough estimate or a lower limit on the size of the broad line region. 
We find that the average mass of the central black holes in blazars is ~ 2.8 x lO 8 Af , with a 
large dispersion, comparable to those computed for other radio-loud active galactic nuclei. 

Key words: galaxies: active — BL Lacertae objects — ultraviolet: galaxies — Gamma rays: 
observations 



1 INTRODUCTION 

Highly Polarized Quasars (HPQ, also referred to as Flat Spectrum 
Radio Quasars) and, occasionally, BL Lacertae objects, collectively 
known as blazars, exhibit broad emission lines superimposed on 
their optical and ultraviolet (UV) continua (Netzer et al. 1994; 
Scarpa, Falomo, & Pian 1995; Vermeulen et al. 1995; Corbett et 
al. 1996; Scarpa & Falomo 1997; Koratkar et al. 1998; Pian et al. 
2002; D'Elia, Padovani, & Landt 2003). Emission lines play an 
important role in the energetics of blazars: some models of multi- 
wavelength blazar emission (Dermer & Schlickeiser 1993; Sikora 
et al. 1994; Ghisellini & Madau 1996) predict that the broad line 
region (BLR) photons are Compton upscattered to X- and gamma- 
ray energies by the relativistic particles composing the jet plasma, 
and form luminous high energy spectral components, which often 
dominate the overall blazar output (Mattox et al. 1997; Bloom et 
al. 1997; Wehrle et al. 1998; Hartman et al. 2001; Ballo et al. 2002; 
Pian et al. 2002). The role of broad line emission in shaping the 
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spectrum of different classes of blazars is however not fully as- 
sessed (Fossati et al. 1998; Ghisellini et al. 1998; Padovani et al. 
2003). 

More in general, the characteristics of the BLR of blazars may 
help in investigating the interplay between the accretion disk and 
the relativistic jet, which is more prominent in blazars than in nor- 
mal QSO and Seyferts (Celotti et al. 1997; Maraschi & Tavecchio 
2003; D'Elia et al. 2003; Wang, Luo, & Ho 2004). AGN broad 
emission lines may be used also to estimate the masses of the 
compact objects residing in the nuclear centers, most likely su- 
permassive black holes (BH), by exploiting their dynamical effect 
on the line-emitting gas clouds. The application of the virial theo- 
rem requires that the size of the BLR is determined either with di- 
rect methods (reverberation mapping technique, Peterson & Wan- 
del 2000; Kaspi et al. 2000; Peterson et al. 2004), or with indi- 
rect arguments (Kaspi et al. 2000; Vestergaard 2002). In the latter 
case the uncertainties are obviously larger, and the methods must 
be tested carefully. Since this has important consequences on the 
evolution and demographics of AGNs, it is crucial to accomplish 
these measurements both for low and high redshift sources. This 
approach has been adopted in the estimate of BH masses of large 
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samples of quasars based on optical emission lines (Woo & Urry 
2002; McLure & Dunlop 2004). The advantage of using UV, rather 
than optical emission lines is that the former correspond to a higher 
ionization state and are therefore presumably more representative 
of the dynamics close to the central massive object. 

In this paper we present the analysis of the broad and intense 
UV emission lines of 16 blazars observed by the Hubble Space 
Telescope (HST) and the Faint Object Spectrograph (FOS). Previ- 
ous studies of AGN UV spectra have been carried out by Bech- 
told et al. (2002), who focussed on the absorption systems, by 
Kuraszkiewicz et al. (2002) and Kuraszkiewicz et al. (2004), who 
present a complete HST FOS atlas of emission line parameters of 
AGNs, and by Evans and Koratkar (2004), who recalibrated the 
pre-COSTAR AGN spectra. We concentrate here on the UV spec- 
tra of blazars. Based on the radiative and kinematic properties of 
the broad emission line region (BLR), we measure the luminosi- 
ties of their BLRs and derive estimates of the BLR sizes and of the 
central BH masses. 

We adopt the "concordant cosmology", fi m = 0.3, Q,a = 0.7, 
and assume Ho = 72 km s _1 Mpc -1 (Spergel et al. 2003). Lumi- 
nosities reported by other authors and used in this paper have been 
transformed into this cosmology. 



2 SAMPLE SELECTION AND DATA ANALYSIS 

We have retrieved from the HST archive 1 all pre- and post- 
COSTAR FOS grating spectra of sources previously classified as 
blazars (Wall & Peacock 1985; Impey & Tapia 1988; Impey & 
Tapia 1990; Impey et al. 1991; Stickel et al. 1991; Stocke et al. 
1991; Padovani & Urry 1992; Wills et al. 1992; Perlman et al. 
1996). We also included in our final list PKS 1229-021, which, 
despite having low polarization (Wills et al. 1992), is considered 
a blazar because of significant emission at MeV-GeV frequencies 
(Hartman et al. 1999), and 3C 273, which has intermediate prop- 
erties between those of blazars (strong radio emission, superlumi- 
nal motion, gamma-ray emission, jet emission dominance at hard 
X-ray energies, e.g., Haardt et al. 1998; Grandi & Palumbo 2004), 
and those of Seyfert galaxies (broad emission lines, big blue bump). 
We selected spectra taken with high-resolution gratings in the UV 
region (G130H, G190H, G270H, G400H). 

This search yielded 24 objects with measurable spectra. Spec- 
tra taken with the same grating within one day were averaged to 
increase the signal-to-noise ratio. For this work we have consid- 
ered only the 16 sources with significant (larger than 3a) emission 
line detections. These are reported in Table 1. 

For 6 sources there are also low-resolution grating (G160L) 
observations in the archive, obtained nearly simultaneously to the 
high-resolution spectra (i.e., within one day). Since the line param- 
eters and spectral indices derived from the former are not signifi- 
cantly different from those measured in the high-resolution spectra, 
we have neglected these spectra. 

Although the considered objects do not represent a complete 
sample, they form a sizeable dataset to investigate the UV line prop- 
erties of blazars. 

After applying a correction for the Galactic absorption using 
the maps of Schlegel, Finkbeiner and Davis (1998) and the extinc- 
tion curve of Cardelli, Clayton and Mathis (1989), we measured 

1 using MAST, the Multi-mission Archive at STScI, see 
http ://archive . stsci . edu 



the equivalent widths (EWs), the intensities and the full width at 
half maximum (FWHM) values of the emission lines fitting a lin- 
ear local continuum on each side of the line (see Table 2). The EW 
uncertainties were estimated by assuming 2a variations of the local 
continuum. 

For each object we have combined the spectra taken quasi- 
simultaneously (within 1 day) with different gratings, excluded the 
regions affected by emission or absorption features, binned the sig- 
nal in 20-50 A wavelength intervals and fitted the continuum with 
a power-law. To account for calibration uncertainties of the data, 
we added a 5% systematic error to the statistical errors. The de- 
rived power-law spectral indices and flux normalizations are given 
in Table 1 . 

A similar analysis of the continuum and line properties of 
these objects has been presented in Kuraszkiewicz et al. (2002) and 
Kuraszkiewicz et al. (2004). 



3 RESULTS 

The redshifts of our objects (see Table 1) range between z = 0.158 
and z = 1.404, with an average value of < z >= 0.84 ± 0.31. 
Thus, the lines typically detected in our FOS spectra are Ly/3, 
Lya, C IV A 1549, C III] A 1909, Si IV A 1400 and in some cases 
Mg II A2798. We report in Table 1 the spectral indices and normal- 
izations of the UV continua and in Table 2 the EW and the FWHM 
values of the emission lines. 

In 3 cases (3C 273, 3C 345, and 3C 454.3) observations at 
more than one epoch are available. Variations of line and contin- 
uum emission are observed, with maximum amplitudes of factors 
of ~2 and 7, respectively. However the observations are too limited 
and sparse to allow a meaningful assessment of correlated line and 
continuum variability. 

In section 3.1 we describe the average properties of the emis- 
sion lines of blazars in the UV spectral region; in section 3.2 we 
use the line and continuum properties of the UV spectra to estimate 
the masses of the central BHs. 



3.1 Emission line properties 

In order to produce a representative high signal-to-noise ratio UV 
spectrum of blazars, we have combined all the UV spectra in our 
dataset. Each spectrum was first reduced to rest-frame and then 
normalized to its average continuum flux. The resulting composite 
blazar spectrum, normalized to unity at the reference wavelength 
of 1500 A, is shown in Fig. 1. The EW, relative intensities, and 
FWHM of the emission lines of the composite spectrum are given 
in Table 3. 

The composite spectrum of blazars is similar to that of nor- 
mal QSO. The line ratios of blazars and normal AGNs are also not 
significantly different. This is illustrated in Fig. 2, where we report 
the luminosities of the Lya and C IV A 1549 lines of blazars, com- 
pared with those of a list of radio-loud quasars (RLQ) observed by 
HST FOS (Wills et al. 1995). Note that in this list of RLQ there 
are 8 objects in common with our sample of blazars; therefore, for 
the purpose of the comparison, these 8 sources have been consid- 
ered as blazars and have been excluded from the RLQ list. We also 
compare with the RLQ 3C 390 (1845+79, z = 0.056), which has 
hybrid properties, i.e., it has substantial polarization (1.3%, Impey 
et al. 1991), but is lobe-dominated at radio wavelengths (Ghisellini 
et al. 1993). Figure 2 shows that the intensity ratio of Lya vs the 
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Table 1. Parameters of Blazars observed with HST FOS 



Object 


Alt. Name 


2 


E B-V 


Date 


Range 6 




«s 






id 

J 1350 


BL0403-1316 




0.571 


0.058 


11.4 Oct 


1991 


1570 


■4780 


1 


20 


± 


0. 


.04 


1.98 


BL0420-0127 




0.915 


0.125 


23.5 Dec 


1996 


2220- 


■3280 





00 


± 


0. 


22 


0.80 


BL0537-4406 




0.896 


0.037 


16.7 Sep 


1993 


2220- 


■3280 





07 


± 


0. 


20 


1.64 


BL0637-7513 




0.656 


0.095 


25.6 May 


1992 


1610- 


■3270 


1. 


25 


± 


0. 


06 


5.64 


BL0954+5537 


4C 55.17 


0.901 


0.0088 


20.9 Jan 


1993 


1570 


■3280 





31 


± 


0. 


08 


0.60 


BL1 144-3755 




1.048 


0.097 


15.6 Jul 


1993 


2220. 


■3280 


0. 


20 


± 


0. 


23 


1.42 


BL1 156+2931 




0.729 


0.019 


26.7 Feb 


1995 


1620- 


■4780 


0. 


52 


± 


0. 


04 


8.29 


BL1226+0219 


3C 273 


0.158 


0.02 


16.8 Jan 


1991 


1090- 


■3280 


1. 


68 


± 


0. 


04 


245 










09.4 Jul 


1991 


1290- 


■3290 


1. 


35 


± 


0. 


06 


117 


BL 1229-0207 




1.045 


0.032 


01.1 Jan 


1995 


1570- 


■3280 


0. 


.81 


± 


0. 


.12 


2.58 


BL1253-0531 


3C 279 


0.538 


0.028 


08.5 Apr 


1992 


1570- 


■4780 


0. 


18 


± 


0. 


03 


1.54 


BL161 1+3420 


DA 406 


1.401 


0.018 


04.9 Apr 


1992 


2220- 


■4780 


0. 


91 


± 


0. 


08 


0.81 


BL1641+3954 


3C 345 


0.595 


0.013 


07.9 Jun 


1992 


1600- 


■4770 


0. 


65 


± 


0. 


04 


3.03 










20.4 Aug 


1995 


1570 


■4800 


0. 


96 


± 


0. 


16 


0.43 


BL2223-0512 


3C446 


1.404 


0.075 


11.8 Sep 


1991 


2220- 


■4780 


0. 


00 


± 


0. 


11 


0.35 


BL2230+1128 


CTA 102 


1.037 


0.072 


12.1 Sep 


1991 


2220- 


■4780 


1. 


13 


± 


0. 


06 


1.95 


BL2243-1222 




0.63 


0.051 


09.4 Oct 


1993 


1610 


■3270 


1. 


59 


± 


0. 


07 


5.18 


BL225 1+1552 


3C 454.3 


0.859 


0.105 


11.9 Sep 


1991 


1620- 


■4770 


0. 


90 


± 


0. 


04 


2.86 










15.3 Nov 


1991 


1620- 


■3270 


0. 


51 


± 


0. 


10 


4.35 










19.4 Aug 


1995 


1700. 


■4800 


1. 


08 


± 


0. 


09 


1.72 



a From the maps of Schlegel et al. (1998). 
b Observed wavelength range, in A. 

c Spectral index of the power-law fitted to the dereddened spectra (fx oc \~ a ). 

d Normalization of the dereddened power-law continuum at 1350 A (rest frame), in 10~ 15 erg s _1 cm~ 2 A -1 . 



Table 2. Emission Line Measurements 



a This table is provided on page 9, in landscape format 



C IV A 1549 in blazars is consistent with that exhibited by normal 
RLQ. 

More in general, we have compared in Table 3 the average 
intensity ratios of the lines we detect in our composite blazar spec- 
trum with those reported by other authors for larger samples of 
QSO or RLQ (Francis et al. 1991; Zheng et al. 1997; Telfer et 
al. 2002). Except for the C III] A1909 line, which appears some- 
what underluminous in blazars, there is good overall agreement. 
This comparison suggests that the structure and physical state of the 
BLR in blazars and normal RLQ are indistinguishable. This is also 
confirmed by the comparison of the Lya line luminosities. In Fig. 3 
we report the continuum luminosity at 1350 A as a function of the 
Lya luminosity for blazars and RLQ. With the exception of 3C 390, 
which is at relatively low redshift, the blazars and RLQ have a sim- 
ilar range of Lya luminosities. The averages of the logarithmic dis- 
tributions, in erg s~\ are < logL(Lya) >= 44.55 ± 0.11 and 
< logL(Lya) >= 44.72 ±0.12 for blazars and RLQ, respectively 
(the uncertainties represent the errors associated with the averages, 
i.e. the standard deviations divided by the square root of the num- 
ber of objects). However, due to relativistic beaming, blazars have 
more luminous continua than RLQ, i.e., blazar lines have smaller 
EW (see Fig. 3). While RLQ emission lines, for any continuum lu- 
minosity, have EW between 100 A and 1000 A, part of the blazars 
exhibit line EW between 10 A and 100 A, and these have the most 



strongly boosted continua. From comparison with the RLQ, we 
have estimated the luminosity enhancement due to beaming. 

Relativistic aberration affects the non-thermal synchrotron lu- 
minosity and depends on the fourth power (for a jet geometry) of 
the relativistic Doppler factor 5 (= [T(l - Pease)]' 1 ). RLQ are 
thought to be the parent population of blazars: their jets are di- 
rected away from the line of sight, so that their luminosities are 
only weakly affected by relativistic beaming. Therefore, we have 
estimated the beaming amplification for the blazars by assuming 
that their continuum luminosities should exhibit a dependence on 
Lya line luminosities similar to that of RLQ (Fig. 3). We fitted the 
RLQ line (in erg s _1 ) and continuum (in erg s _1 A -1 ) luminosities 
to a power-law and obtained the dependence: 

L(1350A) = 0.46 x 10~ 3 L(Lya) 102 , (1) 

which has a scatter of 0.2dex in L(1350 A). We will use this 
relationship in Section 3.2.1 to correct the continuum luminosities 
of blazars for the beaming. 

3.2 Black hole masses 

Under the assumption that the dominant mechanism responsible 
for the width of the broad emission lines is the gravitational poten- 
tial of the central supermassive BH, and that the line widths reflect 
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Table 3. Composite ultraviolet blazar spectrum. 



Line 


EW a 


FWHM a 


Rel. Intensity 6 


Ratio c 


Ratio (LBQS) d 


Ratio (RLQa) e 


Ratio (RLQb)-f 


Ly/3 


6.5 ± 0.5 


10 


8.9 ±0.5 


10.6 


9.3 


19 


19.1 


Lya 


71 ±3 


22 


84 ±3 


100 


100 


100 


100 


SilV 


5.2 ± 0.7 


21 


5.6 ±0.7 


6.6 


19 


6.8 


8.6 


CIV 


45 ± 2 


20 


45 ± 1 


53 


63 


66 


52 


CHI] 


7±1 


23 


6.6 ±0.7 


7.8 


29 


11 


13.2 


Mg II 


19 ±2 


30 


16 ± 1 


19 


34 


24 


22.3 



a At rest frame, in A. Uncertainties are about ~15%. 

b Obtained from EW and continuum normalized to the flux at 1500 A (rest frame). 

c Average percentage intensity ratio with respect to Lya as resulting from our measurements. 

d Same as in Col. 5 for 718 objects in the Large Bright Quasar Survey (Francis et al. 1991). 

e Same as in Col. 5 for for 60 radio loud quasars (Zheng et al. 1997). 

i Same as in Col. 5 for 107 radio loud quasars (Telfer et al. 2002). 
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Figure 1. Composite UV spectrum of blazars obtained from the average of 16 spectra. The most prominent emission lines are labeled. 



the Keplerian velocities of the line-emitting material in a virialized 
system (Wandel et al. 1999; McLure & Dunlop 2001), the BH mass 
Mbh is given by: 



Me 



G 1 v 2 Rblr 



(2) 



where v is the velocity of the gas gravitationally bound to the 
central BH, Rblr is the size of the BLR, and G is the gravitational 
constant. The velocity v can be obtained directly from the FWHM 
of the broad emission lines (v = / x vfwhm), where / is a fac- 



tor that depends on the geometry and kinematics of the BLR (e.g.. 
McLure & Dunlop 2002; Vestergaard 2002). 



3.2.1 Size of the BLR 

The most reliable method to derive .Rblr is through the reverber- 
ation mapping technique (e.g. Peterson et al. 2004, and references 
therein). This uses the time lag of the emission line light curve with 
respect to the continuum light curve to determine the light crossing 
size of the BLR in AGNs. However, this method requires inten- 
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Figure 2. Luminosities of C IV A 1549 and Lyct emission lines for the 
blazars in our list (filled circles). Line measurements of a given source at 
multiple epochs have not been averaged. The Lyct and C IV A 1549 emission 
line intensities have been normalized to their respective average values. The 
solid line is the expected C IV A 1549 to Lyct intensity ratio (L(C IV A 1549) 
= 0.63 xL(Lyct), e.g. Francis et al. 1991). The dashed line is the least square 
regression line L(C IV A1549) = 8.3 xL(Lyct) - 73 ). For comparison, the 
line luminosities of 19 RLQ (open circles) from Wills et al. (1995) and of 
3C 390 (triangle) are also shown. 

sive monitoring of the UV continuum and of the lines and can be 
used only for a limited number of objects (e.g., Korista et al. 1995; 
Onken et al. 2002), including one of our sources, 3C 273 Paltani 
& Tiirler 2005). An alternative way to estimate .Rblr is to use the 
empirical relationship found between Rblr and the optical contin- 
uum luminosity (Kaspi et al. 2000). 

We have derived this relationship in the UV for a sample of 15 
PG quasars and 10 Seyfert 1 galaxies having BLR radii determined 
via reverberation mapping in the optical (Kaspi et al. 2000) and 
measured UV spectra (Vestergaard 2002). One of the PG quasars is 
3C 273, which is also a member of our blazar sample. For this ob- 
ject, we corrected the continuum luminosity for the beaming effect 
(see Eq. 1). For one of the Seyferts, NGC 4151, the UV continuum 
of which varies with high amplitude, we re-measured the contin- 
uum luminosity at 1350 A from the average IUE spectrum. These 
quantities are reported in Fig. 4. We fitted the BLR radii and the 
luminosities at 1350 A (rest frame) with a power-law and obtained 
the following relationship: 

/\r fi A \ °' 61±0 ' 02 
R BLR = (22 .4 ± 0.8) ( ,X g g / ) lt - d ^ S (3) 

If we exclude 3C 273 from the fit, the result is unchanged (the 
power-law index is 0.60 ± 0.02 in this case). We note that the 
slope in Eq. 3 is consistent with that determined in the same wave- 
length range by Kaspi et al. (2005, index 0.56 ± 0.05), it is slightly 
flatter than that found in the optical by Kaspi et al. (2000, index 
0.70 ± 0.03) and slightly steeper than that found by McLure and 



Figure 3. Lyct luminosity vs continuum luminosity at 1 350 A (rest frame) 
for blazars (filled circles) and RLQs (open circles). The triangle represents 
the low redshift RLQ 3C 390. The lines are the loci of the EW = 10 A 
(dashed), 100 A (dotted) and 1000 A (solid). 



Jarvis (2002, index 0.50 ± 0.02) at 3000 A, based on a very similar 
sample of PG quasars and Seyfert galaxies. 

Since the relationship between the BLR radius and continuum 
luminosity is supposed to be valid in the case of a thermal con- 
tinuum (see also discussion in Paltani and Tiirler 2005), we must 
correct the blazar UV continuum luminosities for the effect of rel- 
ativistic beaming. 

For the blazars with continuum luminosity exceeding the 
power-law dependence between the RLQ line and continuum lumi- 
nosities (Eq. 1), we adopted the continuum luminosities computed 
with Eq. 1 at the corresponding line luminosity, and derived the 
BLR radii through Eq. 3. These are reported in Table 4. The cor- 
rection of the continuum luminosity is relevant (i.e., larger than ~3 
times the scatter) for 4 objects (see Fig. 3). We note that our esti- 
mate of the BLR radius of 3C 273 is consistent with that reported 
by Paltani and Tiirler (2005). 

An alternative, independent method for evaluating -Rblr con- 
sists in coupling the luminosity of the BLR with the informa- 
tion carried by the multiwavelength spectrum of the blazar. Since 
blazars, among all AGNs, are the only ones with a spectrum ex- 
tending to gamma-rays, this method is specific for the blazar class 
of AGNs. 

Ten of our blazars have multiwavelength energy distributions 
which have been fitted with synchrotron and inverse Compton ra- 
diation components (Ghisellini et al. 1998). The latter component 
dominates at the X- and gamma-ray energies and originates from 
the scattering of relativistic electrons off both synchrotron pho- 
tons (internal to the jet) and external radiation fields. These include 
broad line photons, the density of which, U ex t, is thus estimated 
through the multiwavelength spectral fit. 

Following the procedure adopted by Celotti, Padovani and 
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(PKS 0537-441, Pian et al. 2002), it appears that this method can- 
not be generalized. 




0.5 - 



42 44 

lc, g( AL i350 l) Terg s->] 



Figure 4. Size of the BLR as a function of the continuum power output 
at 1350 A (rest frame) for the sample of PG quasars (filled circles) and 
Seyfert 1 galaxies (open circles) of Kaspi et al. (2000). The encircled point 
on the top right represents 3C 273. The luminosities have been taken from 
Vestergaard (2002). The line fitting the 2 quantities (Eq. 3) is also reported 
(dotted). 



Ghisellini (1997), we reconstructed the total luminosity of the BLR 
for each of our sources by using the intensities of the observed UV 
emission lines and by assuming for the unobserved lines the line 
ratios of an average quasar spectrum (Francis et al. 1991). These 
derived BLR luminosities are reported in Table 4. 

From the fitted densities of the external photons U ex t and from 
the observed BLR luminosities, the size of the BLR, Rblr, can be 
derived according to: 

where 8 is the relativistic Doppler factor required by the mul- 
tiwavelength modeling. The BLR radii computed with Eq. 4 are re- 
ported in Table 4. We have identified this second method as "spec- 
tral energy distribution (SED) method", in order to distinguish it 
from the one based on the empirical determination of Rblr from 
the continuum luminosity. No clear correlation is found between 
the BLR radii determined with the two methods. One probable ex- 
planation of the discrepancy is that the radiation density of the BLR 
is generally smaller than the parameter U ex t obtained with multi- 
wavelength fits. This parameter includes not only the BLR photons, 
but also additional contributions, such as photons coming directly 
from the accretion disk, or produced by the dusty torus, or by larger 
regions of the jet (Ghisellini priv. comm.). Thus, the SED method 
may underestimate the BLR sizes (and therefore the BH masses) in 
some cases. Moreover, the uncertainties associated with the SED 
fit parameters are large. Therefore, although we had proposed the 
"SED" method for BH mass determination in an individual source 



3.2.2 Mass estimates 

In order to evaluate the BH masses of our objects we have used 
Eqs. 2 and 3, adopting a standard value of / = \/3/2 for the 
kinematic factor, corresponding to an isotropic distribution of the 
BLR clouds (Wandel 1999; Kaspi et al. 2000; Vestergaard 2002), 
that yields virial BH masses consistent with those derived from the 
MsH-Lbuige relationship (Labita, Falomo, & Treves, in prep.). Af- 
ter setting vfwhm to suitable units we obtain the relation: 

3.26x10 ^ 10 44 ergs -i ) 1 103 kms-ij Mq(5) 

Vestergaard (2002) also derives a formula for the central BH 
mass, based on the continuum measurements at the rest frame 
wavelength of 1350 A of a sample of 26 AGN with BLR radii de- 
termined by reverberation mapping. However, this relationship was 
calibrated against the BH masses determined from optical measure- 
ments, and not directly from the BLR size, as we do. 

By using Eq. 5 with the beaming-corrected L(1350 A) lumi- 
nosities and vfwhm estimated from our spectra, we have computed 
the central BH masses. For consistency with Vestergaard (2002) we 
have used the FWHM of C IV A1549. These masses are reported 
in Column 5 of Table 4, and in Fig. 5, where they are compared 
with those computed with our Eq. 5 for a sample of PG quasars, for 
which Vestergaard (2002) reports UV luminosities and C IV A 1549 
line FWHM values. The blazar BH masses are statistically con- 
sistent with those of quasars: the averages of the logarithmic dis- 
tributions, in solar masses, are < logMeH >= 8.31 ± 0.10 and 
< logMeH >= 8.42±0.08 forblazars and RLQ, respectively (the 
uncertainties represent the errors associated with the averages, i.e. 
the standard deviations divided by the square root of the number of 
objects). 

For comparison, we have also computed the blazar BH masses 
using Vestergaard's relationship (Eq. 8 of Vestergaard 2002), 
with our measured blazar luminosities and FWHM values of the 
C IV A 1549 line. These are systematically lower than the corre- 
sponding ones determined using Eq. 5, by a factor 1.4-2. We have 
also compared our estimated BH masses with those determined for 
the same blazars by other authors: a number of our objects are in 
common with the samples of Liang and Liu (2003), Woo and Urry 
(2002), and Wang et al. (2004). Their mass estimates have been re- 
ported in the last 3 columns of Table 4: our masses are generally 
smaller. This may be partially due to our correction of the contin- 
uum luminosity for the beaming effect. In particular, this must be 
the case for two of the four objects in our sample with the most 
strongly beamed continuum, 3C 279 and 0537-441. 

We do not find a correlation of the BH mass with redshift for 
our blazar sample. 



4 SUMMARY AND CONCLUSIONS 

We have studied the properties of the UV emission lines of blazars, 
mostly from single epoch HST FOS spectra, and found that the av- 
erage blazar UV spectrum is similar to that of RLQ. This is the sum 
of a thermal and non-thermal component. Our targets are mainly 
HPQ and Low-Frequency Peaked BL Lacs (Padovani & Giommi 
1995; Fossati et al. 1998) where the emission of the non-thermal 
synchrotron component peaks at optical/IR frequencies. Therefore 
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Table 4. BLR Luminosities and Central Black Hole Masses. 



UDject 


r a 
BLR 


n BLR 


nc 

BLR, SED 


71 /f d 
BH 


BH ,LL 


A/T f 

BH,WU 


BH,W 


BL0403-1316 


22.6 


145 




2.4 




12 




BL0420-0127 


18.7 


145 


438 


2.3 


8 


11 


9 


BL0537-4406 


6.93 


77 


211 


0.5 


16 




5 


BL0637-7513 


50.2 


285 




1.9 




26 




BL0954+5537 


5.38 


66 


312 


0.5 


8 


1.2 




BL1 144-3755 


2.45 


34 




0.4 








BL1 156+2931 


13.7 


115 


158 


4.3 


8 






BL1226+0219 


33.8 


219 


17 


4.0 


0.2 




16 


BL1229-0207 


74.1 


357 


117 


7.4 


10 






BL1253-0531 


2.42 


36 


104 


0.8 


4 


2.7 


3 


BL161 1+3420 


34.1 


196 


92 


6.2 


40 


37 




BL1641+3954 


12.5 


92 




2.7 




26 


28 


BL2223-0512 


25.2 


150 




2.8 






6 


BL2230+1128 


41.4 


260 


33 


3.1 


12.6 






BL2243-1222 


47.8 


271 




2.1 








BL2251+1552 


33.3 


223 


29 


3.1 


25 


15 


13 



a Luminosity of the BLR, in units of 10 erg s . 
6 Size of the BLR computed via Eq. 3, in light days. 
c Size of the BLR computed using the SED method, (Eq. 4), in light days. 

d Mass of the central BH computed using Eq. 5, in units of 10 8 Mq. The statistical uncertainties, dominated 

by the scatter of Eq. 3, are about a factor 2. 
e Mass of the central BH from Liang & Liu (2003). 
f Mass of the central BH from Woo & Urry (2002). 
f Mass of the central BH from Wang et al. (2004). 
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Figure 5. Histograms of central BH masses of AGNs. The simply hatched 
area represents the masses of PG quasars computed with Eq. 5 from the 
luminosities and C IV A 1549 emission line FWHM values reported by from 
Vestergaard (2002). The double-hatched area represents the BH masses of 
our blazar sample. Mass estimates obtained for a same source at different 
epochs have been averaged. 



a large relative contribution from the thermal accretion disk is ex- 
pected in the UV region (e.g., Bregman et al. 1986). This is clearly 
evident from the spectral energy distribution of 3C 273 (Ulrich et 
al. 1980; Courvoisier 1998) and may be significant in 3C 279 (Pian 
et al. 1999). 

With the aim of estimating central BH masses of blazars, we 
have assumed Keplerian conditions in the BLR gas motion and 
have evaluated the BLR size using the results of a fit of UV lu- 
minosities and BLR radii of a sample of QSOs having BLR sizes 
determined via reverberation mapping in the optical. We have de- 
rived a relationship between Rblr and luminosity in the UV do- 
main (1350 A at rest frame), which exhibits a slope consistent with 
that of Kaspi et al. (2005), although slightly flatter and steeper than 
proposed in the optical and near-UV by Kaspi et al. (2000) and 
McLure & Jarvis (2002), respectively. 

For those 10 blazars having multi wavelength spectral fits we 
have also applied an independent method of BLR size determina- 
tion, based on the combination of the observed Lblr and fitted ex- 
ternal radiation density ("SED" method). We have not found a clear 
correlation between the BLR sizes obtained with the two meth- 
ods, with the largest deviations observed in the sense of a deficit 
of the "SED" radii with respect to those obtained with the empiri- 
cal Rblr — AL^ relationship. We conclude that the SED method 
yields a BLR size inconsistent with that derived from the contin- 
uum luminosity. 

Our estimated BH masses have an average of (2.8 ± 2.0) x 
10 Mq (the quoted uncertainty is the standard deviation) and are 
comparable with those of lower redshift blazars, estimated with dif- 
ferent methods (Barth, Ho, Sargent 2003; Falomo, Carangelo, & 
Treves 2003). 

The distribution of our blazar BH masses computed with Eq. 5 
is consistent with the distribution of the PG quasar masses com- 
puted with the same equation. These results suggest that the differ- 
ences between radio powerful sources and radio-weak ones are not 
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due to the mass of the BHs residing at their centers. However, the 
validity of this conclusion at the intermediate/high redshifts must 
be corroborated by the analysis of wider samples of homogeneous 
datasets. Moreover, further intensive spectroscopic monitoring of 
the brightest blazars at optical and UV wavelengths is required, in 
order to construct well sampled continuum and emission line light 
curves for the application of the reverberation mapping technique. 
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Table 2: Emission Lines Measurements 



Object 


z 


Ly/3° 




Lya" 




Si 


IV 


C IV 




C III 




Me 


II 






EW 6 FWHM 6 


EW 


FWHM 


EW 


FWHM 


EW 


FWHM 


EW 


FWHM 


EW 


FWHM 


BL0403-1316 


0.571 






135 ±4 


14 


9.8 ±0.7 


22 


128 ±4 


17 


26 ± 1 


18 


33 ±3 


30 


BL0420-0127 


0.915 






105 ±6 


11 


7.2 ±0.8 


19 


56 ±2 


17 










BL0537-4406 
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20 ±2 
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1.7 ±0.4 


13 


11.4 ±0.7 


11 
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11 ±2 
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89 ±2 


11 
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43 ±2 
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12 ±1 


13 
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40 ± 1 


ll c 






25 ±1 
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4.6 ±0.8 


14 










BL1156+2931 
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12.0 ±0.8 
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1.1 ±0.3 
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8.5 ±0.5 


26 
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0.158 


2.6 ±0.5 
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52 ±2 


12 






26.8 ±0.9 


18 


8.8 ±0.3 


19 


6.6 ±0.5 


24 
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0.158 






43 ±3 
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31 ±2 
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9.1 ±0.6 
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8.5 ±0.6 


25 


BL1229-0207 


1.045 


12 ±1 


13 


94 ±4 


13 


5.7 ±0.5 


25 


42 ±2 
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BL1253-0531 
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22 ±1 
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20 ± 1 


20 c 


3.9 ±0.5 
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6.4 ±0.5 
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14.8 ±0.9 


21 


22 ± 1 
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149 ± 36 
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16 ±5 
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150 ± 16 
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45 ±8 
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91 ±19 
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1.404 


23 ±2 


15 


120 ±3 
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58 ±5 


19 


14 ±3 


40 






BL2230+1128 


1.037 






77 ±2 


13 


2.9 ±0.4 


19 


29 ± 1 
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21 
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13 
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15 ±1 
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89 ±2 
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3.2 ±0.5 
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31.8 ±0.8 


16 


13.4 ±0.7 


24 
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" The Lya and Ly/3 line EW measurements may have a residual contamination by the N V A1240 and O VI A1034 lines, respectively. This possible 

contamination does not exceed ~10%. We have deblended the contamination in the FWHM measurements. 
6 At rest frame, in A. Uncertainties on the FWHM values are about 15%. 
c This line has an asymmetric profile, with a more prominent red wing. 
d 16.8 Jan 1991. 
e 9.4 Jul 1991. 
f 7.9 Jun 1992. 
i 20.4 Aug 1995. 
h 11.9 Sep 1991. 
i 15.3 Nov 1991. 
i 19.4 Aug 1995. 



